Inflation in the warm and cold regimes 
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It is now understood that inflation dynamics conies in two forms, isentropic or cold inflation 
and nonisentropic or warm inflation. In the former, inflation occurs without radiation production, 
whereas in the latter both radiation production and inflation occur concurrently. Recent, detailed, 
quantum fleld theory calculations have shown that many generic inflation models, including hybrid 
inflation, which were believed only to have cold inflation regimes, in fact have regimes of both 
warm and cold inflation. These results dispel many foregone assumptions generally made up to now 
about inflation models and bring to the fore various elementary issues that must be addressed to do 
reliable calculations from inflation models. Here I review these results and issues. I then show that 
warm inflation has intrinsic model independent features that makes it natural or equivalently have 
no "eta problem". Next density perturbations and observational consequences of warm inflation 
are discussed. Finally the implications of warm inflation to model building and physics beyond the 
Standard Model are outlined. 

Plenary talk at Cosmion04. appears in conference proceedings. 
Gravitation and Cosmology, V. 11, No. 1-2 (41-42), 51 (2005). 



I. INTRODUCTION 

For over two decades, inflation has been a very success- 
ful idea. In its earliest days, this success was attributed 
to the ability for this idea to unite particle physics and 
cosmology. In the past decade, the success of the in- 
flation idea has been driven by its consistency to obser- 
vation based on data in particular from precision CMB 
satellite experiments. The growing success over the years 
of the inflation idea has also led to an increasing under- 
standing of the underlying dynamics of inflation. By now 
it is known that there are two dynamical realizations of 
inflation. One is the original or standard picture, also 
referred to as isentropic or cold inflation [E H 01 ■ In this 
picture inflationary expansion occurs with the universe 
in a supercooled phase which subsequently ends with a 
reheating period that introduces radiation into the uni- 
verse. The fluctuations created during inflation are effec- 
tively zero-point ground state fluctuations and the evo- 
lution of the inflaton field is governed by a ground state 
evolution equation. In an alternative class of cold infla- 
tion models, inflation has a geometrical origin 4] with 
adiabatic density perturbations Q- The other picture 
of inflation dynamics is nonisentropic or warm inflation 
p]. In this picture, inflationary expansion and radiation 
production occur concurrently. Moreover, the fluctua- 
tions created during inflation emerge from some excited 
statistical state and the evolution of the inflaton has dis- 
sipative terms arising from the interaction of the inflaton 
with other flelds. 

The dividing point between warm and cold inflation is 
roughly at pr ~ H, where pr is the radiation energy 
density present during inflation and H is the Hubble pa- 
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rameter. Thus pr > H is the warm inflation regime 
and Pr < H is the cold inflation regime. This criteria 
is independent of thermalization, but if such were to oc- 
cur, one sees this criteria basically amounts to the warm 
inflation regime corresponding to when T > H. This is 
easy to understand since the typical inflaton mass dur- 
ing inflation is m^ « H and so when T > H, thermal 
fluctuations of the inflaton field will become important. 
This criteria for entering the warm inflation regime turns 
out to require the dissipation of a very tiny fraction of 
the inflaton vacuum energy during inflation. For exam- 
ple, for inflation with vacuum {i.e. potential) energy at 
the GUT scale ~ 10^^~^^GeV, in order to produce ra- 
diation at the scale of the Hubble parameter, which is 
w 10^"~^"'^GeV, it just requires dissipating one part in 
10^" of this vacuum energy density into radiation. Thus 
energetically not a very significant amount of radiation 
production is required to move into the warm inflation 
regime. In fact the levels are so small, and their eventual 
effects on density perturbations and inflaton evolution 
are so significant, that care must be taken to account for 
these effects in the analysis of any inflation models. 

In this paper, we will examine warm inflation dynamics 
and compare it to cold inflation dynamics. In Sect. HTIthe 
basic equations of both pictures are reviewed. In Sect. 
mil a particular mechansim for dissipation is presented 
and the dissipative quantum field theory calculation is 
done to determine dissipative effects in the inflaton effec- 
tive equation of motion. In Sect. IIVI densitv fluctuations 
in cold and warm inflation are presented. A calculation 
of inflation for a momomial potential is done, and it is 
shown that in contrast to cold inflation, in the warm in- 
flation case the mass of the inflaton is bigger than the 
Hubble parameter and the field amplitude is below the 
Planck scale. These are two very different outcomes of 
warm infiation that are not found in the cold infiation 
case for such potentials. In Sect. Elan infiation analysis 
of the SUSY hybrid model is done. The parameter space 



is shown to divide into regimes of cold and warm infla- 
tion. The spectral index is also calculated to examine 
observational differences in the two cases. Finally some 
examination is made of particle physics of this model in 
the warm inflation case. Finally in Sect. IVf I our conclu- 
sions are given. 



II. TWO INFLATION DYNAMICS 

There are two distinct dynamical realizations of infla- 
tion. In the original picture, termed cold, supercooled or 
isentropic inflation Q, |^ y| , the universe rapidly super- 
cools during inflation and subsequently a reheating phase 
is invoked to end inflation and put the universe back into 
a radition dominated regime. In this picture, inflaton 
dynamics is governed by the equation 



(p + ilKp + CR(^ - 
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where R = 6{d/a + a? /a^) is the curvature scalar. The 
conditions for slow-roll inflation in this picture require 
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where nip =— ■mpi/\/^ — 2.4 x 10^^ GeV is the reduced 
Planck mass. In this slow-roll regime the evolution equa- 
tion (^ is well approximated by 



dp 



(4) 



In the other picture, termed warm or nonisentropic in- 
flation JJJ , dissipative effects are important during the in- 
flation period, so that radition production occurs concur- 
rently with inflationary expansion. Phenomenologically, 
the inflaton evolution in simple warm inflation models 
has the form, 



^+[2,H + T{ip)]p + i'RLp- 
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For T = 0, this equation reduces the the familiar inflaton 
evolution equation for cold inflation Eq. J^l, but for a 
nonzero T, it corresponds to the case where the inflaton 
field is dissipating energy into the universe, thus creating 
a radiation component. The conditions for slow-roll in- 
flation are modified in the presence of the extra friction 
term T, and we have now: 
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where r = T/{3H), and en, rjn are the slow-roll param- 
eters without dissipation given in Eqs. (0) and JSJ. In 
addition, when the friction term T depends on the value 
of the inflaton field, we can define a third slow-roll pa- 
rameter 
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In this slow-roll regime of warm inflation the infl;aton 
evolution equation is well approximated by 



[3H + Tiip)]p + ^Rp- 
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The dissipation of the inflaton's motion is associated 
with the production of entropy. The entropy density of 
the radiation s{4>^T) is defined by a thermodynamic re- 
lation in terms of the thermodynamic potential. 
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The rate of entropy production can be deduced from the 
conservation of energy-momentum. The total density p 
and pressure p are given by. 
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Energy-momentum conservation, 

p + iH{p + p)^{), 



(12) 
(13) 

(14) 



now implies entropy production. Making use of Eq. © 
(with ^ = 0) we get 



T{h + iHs) = T(/)' 
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(15) 



The zero curvature Friedman equation completes the set 
of differential equations for 0, T and the scale factor a. 



3iy2 = 87rG(i(/)2 + F + Ts). 



(16) 



The entropy production has been described in a 
slightly different way in the initial warm inflation papers 
ja, 0] ■ We can recover an alternative equation in the case 
when the temperature corrections to the potential are 
negligeable. If we set 5mT = in the finite temperature 
effective potential 14// {(j), T) , then the radiation density 
Pr = 4:sT/3, and Eq. H15|l becomes 



i2 



iHpr =T<j)\ 
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This equation is only valid when dniT = 0. 

Eq. H1U|) is the simplest form of warm inflation dy- 
namics. The basic idea of warm inflation is that radiation 



production is occurring concurrently with inflationary ex- 
pansion due to dissipation from the inflaton field system. 
This dissipation would imply nonconservative terms in 
the inflaton evolution equation. In gneeral these terms 
need not be of the simple temporally local form Tip as 
in Eq. HlOfl . but could be nonlocal. One such example 
would be a form like J K{t, t')if"-{t')dt' . Moreover in gen- 
eral this radiation production could be produced far from 
equilibrium. 



and everything else is the environment, which in partic- 
ular means the $ fluctuation modes, the scalars Xj ^-nd 
the spinors il^k , ipk are regarded as the environment bath. 
In a Minkowski background, at T = 0, the EOM for (p 
has been derived in ^ using the Schwinger closed time 
path formalism. Here we follow a completely analogous 
approach and derive the EOM in a FRW background 
following [31 . The field equation for $ can be readily 
obtained from Eq. (|19|l and it is given by 



III. RADIATION PRODUCTION MECHANISM 
DURING INFLATION 

The key question is what types of first principles mech- 
anisms can result in radiation production during infia- 
tion. Here we consider the example of [3, 13, which devel- 
ops the mechanism of the scalar inflaton field (j) exciting a 
heavy bosonic field x which then decays to light fermions 
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This mechanism is expressed in its simplest form by an 
interaction Lagrangian density for the coupling of the 
inflaton field to the other fields of the form 



>C/ = -7;9^^^X^ 



g'^^x^x - hx^d^d, (19) 



where f/'d are the light fermions to which x-particles can 
decay, with m^ > 2m^^ . Aside from the last term in Eq. 
(|19|l . these are the typical interactions commonly used 
in studies of reheating after inflation [l^, U^ ■ However a 
realistic inflation model often can also have additional in- 
teractions outside the inflaton sector, with the inclusion 
of the light fermions V'd as depicted above being a viable 
option; an example of this is the SUSY hybrid model for 
which a recent study of warm inflation has been done [IJ 
and will be discussed in Sect. Moreover in minimal 
supersymmctry (SUSY) extensions of the typical reheat- 
ing model or multifield infiation models, the interactions 
of the form as given in Eq. (|19(l can emerge as an au- 
tomatic consequence of the supersymmetric structure of 
the model. Since in the moderate to strong perturbative 
regime, reheating and multifield inflation models will re- 
quire SUSY for controlling radiative corrections, Eq. H19|) 
with inclusion of the ipd field thus is a toy model repre- 
sentative of many realistic situations. 

The effective evolution equation for the inflaton back- 
ground field arising from mechanism Eq. I|19|l has been 
computed in 8, 9J. The basic idea is we are interested in 
obtaining the effective equation of motion (EOM) for a 
scalar field configuration ip = ($) after integrating out 
the $ fluctuations, the scalars Xj and spinors •ipk,'>Pk- 
This is a typical " system-environment" decomposition of 
the problem in which ip is regarded as the system field 
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In order to obtain the effective EOM for (p, we use the 
tadpole method. In this method we split 4> in Eq. if^ . 
as usual, into the (homogeneous) classical expectation 
value ip{t) = ($) and a quantum fluctuation 4>{x), $(a;) = 
ip{t) + 4>{x). In this case, the field equation for $, after 
taking the average (with {(j>{x)) — 0), becomes 



m + 3^ipit)+mlipit) + ^ip^t) 

+ ^i?(t)^(t) + ^^(t)(02) + ^(03^ 

+ Efff bW(4) + ('^X,')]=0, (21) 

where ((/>^), {(f)^), (x?) and ((/ix^) can be expressed J8| in 
terms of the coincidence limit of the (causal) two-point 
Green's functions Gt^{x,x') and G^^{x,x'), for the <I> 
and Xj fields respectively. In particular 

G++{x,x')=i{T+Xjix)Xj{^')) 
Gx^ix,x) = i(T^Xj ix)Xjix')) 
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and similarly for G^{x,x'). The momentum space 
Fourier transfor of G^ (x, x') is 



G^. [x,x ) — i 
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where 
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Gj+(q,t,i') = G>.(q,t,i')e(t-t') 
+ G<.(q,t,i')e(t'-i), 
G--(q,t,0 = G>.(q,t,t')e(i'-i) 



+ G<.(q,t,t')^(i- 
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G+-{q,t,t') = G< (q,t,0, 




G-+(q,t,i') = G>.(q,t,i') 





(24) 

and similarily for G^{x,x'). The Green's functions 
G?':^(q, t, i') are written in terms of the modes of the 
scalar field as 

G>(q,i,<') = /^^,i(q,t)/^^;2(q,t')^(i-i') 

G<(q,t,i') = /*^^i(q,t)/*^^2(q,t')e(t-i') 

+ /x,;i(q,t')/x,;2(q,i)^(^'-^), (25) 

where in general there will be two independent solutions 
/x3;i,2(q, i) since Xj obeys an equation second order in 
time. To obtain the evolution equations for the mode 
functions in Eq. [23 the fermion fields that interact with 
Xj fields are integrated out following Q , which then leads 
to 






,d d 
a dt 
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+ j di'a3(i')n^, (q; i, i')/x. (q. ^') - , (26) 

where li^ . (q; i, t') is the spatial Fourier transform of the 
Xi field self-energy term. In particular, in terms of the 
self-energy matrix on the closed time path 11^^^ is defined 
through 



I]++(x,x') + S+-(x,x') ^^ 
+- 

Xi 
Xj 



Xj 
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ni,x^(a;,a:') = [20(ti - is) - 1] 
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which have the properties Hij^, (^^j 2:') = IIi.^^. (x', x) and 

Typically, equations for the mode functions for an in- 
teracting model, of the general form as given by Eq. (|26|l , 
can be very difficult to solve analytically, in particular 
for an expanding background. There are a few particu- 
lar cases, such as for de Sitter expansion H ~ constant, 
so a{t) = exp{Ht), and power law expansion a{t) ~ t", 
where solutions for the mode equation for free fluctua- 
tions are known in exact analytical form (see e.g. Ref. 
|13J|). However, for deriving an approximate solution for 
the mode functions in the interacting case, we can apply 



a WKB approximation for equations of the general form 
Eq. (|26|l and then check the validity of the approxima- 
tion for the parameter and dynamical regime of interest 
to us. As will be seen below, under the dynamical con- 
ditions we are interested in studying in this paper, this 
approximation will suit our purposes. Let us briefly recall 
the WKB approximation and its general validity regime, 
when applied to obtaining approximate solutions for field 
mode equations. An approximated WKB solution for a 
mode equation like 



r w2 



dt^ 



+ ^'(q,i) 



/(q,t)-0 



(29) 



is of the form 
fwKB{(i,t) = l/[w(q,i)]^/^exp 



±i 



dt"uj{q,t") 



(30) 

which holds under the general adiabatic condition 
w(q,t)<tj2(q,t). 

Proceeding with our derivation, consider then a dif- 
ferential equation in the form of Eq. H26|) . Instead of 
working in cosmic time, it is more convenient to work in 
conformal time r, defined by dr — dt/a{t), in which case 
the metric becomes conformally fiat. 



ds2 ^ a(^)2 ^^^2 _ ^^2 



(31) 



By also defining a rescaled mode field in conformal time 

by 
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/(q,T)=/(q,t), 



(32) 



we can then re-express Eq. (|26ll in the form (generically 
valid for either (f> or Xj scalar fiuctuations) 



-^/(q,r)+cj(q,T)2/(q,^) 
+ ydT'n(q,r,r')/(q,r') = 0, 

where we have defined 
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In H34|l the conformal symmetry appears in an explicit 
form, with S, — 1/6 referring to fields conformally cou- 
pled to the curvature, while ^ = gives the minimally 
coupled case. Note also that in conformal time the scalar 
curvature becomes 



R{t) 



6 d^g 



(35) 



In Eq. I|33|l we have also defined the self-energy in con- 
formal time as, 



which is then used in Eq. H40|l for the second derivative 
term to determine the next order approximation, 



n(q, 



a(T)3/2a(T')3/2 



= n(q,i,t') 



(36) 



where the self-energy contribution H, coming from the 
integration over the bath fields, is given by the space 
Fourier transformed form for Eq. H27|l . In H27|l . 11 was 
split into symmetric and antisymmetric pieces with re- 
spect to its argument as defined in Eq. H28|l . Thus based 
on Eq. (|28l) , the self-energy term in H33|l can then be writ- 
ten as n(q,T,T') = ni(q, r, r')-|-n2(q, r, r'). In addition, 
by writing the self-energy term in a diagonal (local) form 



n(q, r,r') = n(q, t)(5(t - r') 

= [ni(q,T)+n2(q,r)](5(T-r'), (37) 

and from the properties satisfied by Hi and 112 , it results 
that ni(q, r) must be real, while Il2(q, t) must be purely 
imaginary. The real part of the self-energy contributes 
to both mass and wave function renormalization terms 
that can be taken into account by a proper redefinition 
of both the field and mass M. On the other hand, the 
imaginary term of the self-energy is associated with de- 
caying processes, as discussed previously. So, we can now 
relate the decay width in terms of the CTP self-energy 
terms as 



imn 



2<I> 2u! 

and Eq. (|33|l can be put in the form 



-^ + (i(q, r)2 - 2iw(q, T)f (q, r) 



(38) 



We now proceed to obtain a standard WKB solution 
for Eq. 1(23 . To do this, following the usual WKB 
procedure, we assume the solution to have the form 
/(q, r) — cexp[i7(q, r)], where c is some constant that 
can be fixed by the initial conditions, given by (|44() be- 
low. This form of the solution is then substituted into 
(|S^ to give 



ij" - 7'^ + cj2 - 2i(I>r = . 



(40) 



Working in the standard WKB approximation, for the 
zeroth order approximation we neglect the second deriva- 
tive term in H40|) . Then, by taking F <C w, we obtain 



70 w T / At' (w - %Y) 

Ira 



(41) 
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dr' [w-if -|-0((I>'V<^^)] +ilnVS. (42) 



The next and following orders in the approximation bring 
higher powers and derivatives of (D'/(D^, which in the adi- 
abatic regime, oj'/a)^ <C 1, are negligible and we are then 
led to the result 



/i,2(q,T) w ^=exp 



TJ / dr' (w - iV) 

Jrn 



(43) 



The solutions for the modes of the form Eq. (|43|l and 
their complex conjugate are general within the adiabatic, 
or WKB, approximation regime of dynamics. Finally, we 
completely and uniquely determine the modes by fixing 
the initial conditions at some initial reference time tq, 
which can be chosen such that in the limit of fc ^ 00 
or iJ — > we reproduce the Minkowski results. These 
conditions, which correspond to the ones for the Bunch- 
Davis vacuum 'iSll, can be written as 



/i,2(q,To) 



V2^^M ' 



/ i,2(q7To) = T«\Mn))72 



(44) 



which already fixes the constant c in (|43|l as c = 1/V2. 

Using the above results in (|25|l and after returning to 
cosmic time i, we obtain the result, valid within the WKB 
approximation, or adiabatic regime. 



G>«)(q,t,0 
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\a{t)a{t')fl^ 



G>«)(q,t,0, (45) 



/(q,T)-0. (39) where 
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-ij^, dt"luj{t")-ir{t")]g,^ _ ^,^ 



G>(q,t,t') 



+ e-^/>*"["^*")+^^^*"^l0(^'-t)| , 
G<(q,i,t') = G>(q,t',i), (46) 

where F is the field decay width in cosmic time, obtained 
from (|38|l and 
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KP{t) 



with Af^(t), for $ particles, given by 



(47) 



Mlit) = m| + ^^{tf + (e - Rit) , (48) 

while for Xj particles, 
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where 



Chiith[f{t),<l>{x),Xj,-ii^k,ipk,gf,,y] , (54) 



K(') 



mz 



-g^tf 



e 



Rit) 



(49) 



The same result Eq. ()46|l could in principle be inferred 
in an alternative way by expressing the propagator ex- 
pressions in terms of a spectral function, defined by a 
Fourier transform for the difference between the retarded 
and advanced dressed propagators, 

G'°\x,x') = e{t-t')[G>[x,x')^G<{x,x')\ 

= G++{x,x')~G+-{x,x') , (50) 

G'^'^^ix^x') = e{t' -t)[G<{x,x')-G>{x,x')\ 

= G++{x,x')~G-+{x,x') , (51) 

and approximating the spectral function as a standard 
Breit-Wigner form with width given by T and poles de- 
termining the arguments of the exponential in (|43|l and 
its complex conjugate J16j. The validity of this approxi- 
mation in particular was recently numerically tested and 
verified in Ref . jl3| for a 1-1- 1 d scalar field in Minkowski 
space-time. In the Minkowski space-time case, results 
analogous to Eg l|^ were explicitly derived in Refs. 
Hi HE HM H^ I2OI ■ Indeed, for the case of no expansion 
a{t) = constant, Eq. (|46|l reproduce the same expressions 
as found in the case of Minkowski space-time. 

The result Eq. (|46|l , from the previous approximations 
used to derive the WKB solution Eq. (|43|l . is valid under 
the requirements 



r^ < W0 , 

and the adiabatic conditions. 
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(53) 



where in the second term in the equations (|53|) we have 
made the change back to comoving time and used Co = 
a{t)io(t) = a^Jq^/a? + KP. 

We now turn our attention to the EOM Eq. (|^ . 
where we will work it out in the response theory ap- 
proximation similar to the treatment in |3| • Consider the 
Lagrangian density in terms of the background (system) 
field tf{t) and the fluctuation (bath) fields. 






i^[tf^^v{tf-\R^{tf 



is the sector of the Lagrangian independent of the fluc- 
tuation bath fields, while £bath denotes the sector of the 
Lagrangian that depends on the bath fields and in par- 
ticular includes the key interaction terms Eq. (|19|1 . In 
the following derivation it will be assumed that the back- 
ground field (p{t) is slowly varying, something that must 
be checked for self-consistency. Thus, if we consider the 
decomposition of ^p{t) around some arbitrary time to as 
(p(t) = (^(to) + Sip{t), Sip{t) can be regarded as a pertur- 
bation, for which a response theory approximation can 
be used for the derivation of the field averages in Eq. 

In response theory we express the change in the 
expectation value of some operator 0{t), 5{0{t)) — 
(C'(i))pert — {0{t)), under the influence of some exter- 
nal perturbation described by Ifpert which is turned on 
at some time ioj as (for an introductory account of re- 
sponse theory, see for instance Ref. |2l] |) 



5m)) 



dt'i 



Hp,,t{t'),d{t) 



(55) 



where the expectation value on the RHS of Eq. (|55|) is 
evaluated in the unperturbed ensemble. The response 
function defined by Eq. (|55|l can be readily generalized 
for the derivation of the field averages. Provided that the 
amplitude S(p{t) is small relative to the background field 
f{to), perturbation theory through the response function 
can be used to deduce the expectation values of the fields 
that enter in the EOM Eq. (|21|l . In this case the perturb- 



Sip 



t6ip 



ing Hamiltonian -ffpcrt is obtained from /^jj^^, where C-^ 
is the part of the interaction Largangian proportional to 
5ip. From £jj^ and Eq. l|^ we can then determine the 
averages of the bath fields, for example (Xj(i))) as an 
expansion in 5Lp{t), starting from the time to and in an 
one-loop approximation, as 
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where 



([X,'(x,i),xf(x,i')]> = 2*Ini(Txf(x,i)x|(x,i')> 



4i 



d3 



Similarly the expression for the other expectation values 
((^^), {(f>^), and (0Xj) can be determined. 

Substituting these field averages into Eq. l(^ we then 
obtain the (^-effective equation of motion 



ip{t) + 3H{t)ip{t) 



dV;^{^{t),R{t)) 
dLp{t) 



+AV(i) / dt' ipit')ip{t')K4t,t') 

Jto 



i=i 



^4gXi) / dt'^it')^it')K^^{t,t') = , (58) 



where V^g stand for the renormalized effective potential, 



KxAt,t') 
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druJxjAT) 
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'i^x.At)^x.An 



(59) 



and similar expression for K^{t, t'). 

In the case where the motion of ^ is slow, the above 
equation can be further simplified. In particular, in this 
case the adiabatic-Markovian approximation can be ap- 
plied to the nonlocal kernels K^.{t,t') etc... This con- 
verts Eq. (|58|l into one that is completely local in time, 
albeit with time derivative terms. The details of this 
approximation for Minkowski space-time can be found 
in 13 . Its extension to an expanding FRW background 
follows analogous lines. The Markovian approximation 
amounts to substituting i' — > i in the arguments of the 
93-fields in the kernels in Eq. (|58|1 . The adiabatic approx- 
imation then requires self-consistently that all macro- 
scopic motion is slow on the scale of microscopic mo- 
tion, thus Lpjip^H < Ty^. Moreover when H < M^, the 
kernel K^{t,t') is well approximated by the nonexpand- 
ing limit H ^ 0. The validity of all these approxima- 
tions were examined in [9j . The result of these approx- 
imations is that, after trivially integrating over the mo- 
mentum integral in the last term in Eq. (|58|l . the effec- 
tive EOM Eq. l(5S|) becomes Eq. Q. By setting the 
couplings gj — 5' = 5, ft-fej = h ^ g ^ X, the mass 
M^ ~ g(^ > m.0^ and F^ ~ N^h'^M'^/lS-Kto^], it leads to 
the friction coefhcient T{(p) in Eq. 



m 



T(<^) = N^ 



\f2g'^a^ip^ 



647rMx^/l-|-a2 



(60) 
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exact 
nonexpanding kernel 

adiabatic-Markovian 
no dissipation 
g = 0.37 
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O 

~o 
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20 
tx[10"'^GeV"^ = 0.45/H(0)] 

FIG. 1: Evolution of if>{t) for A = 10-l^ g = h = 0.37, C = 0, 
V5(0) = mpi, (^(0) = 0. 



Fig. ^compares the various approximations for a rep- 
resentative case where g = h = 0.37 and the inflaton 
potential is that for chaotic inflation Vcs{(p) = Xip^/A 
with A = lO""'^^ I3]. In Fig ^ evolution has been exam- 
ined at the final stages of chaotic inflation where we start 
with ip{to = 0) = mpi. The solid line is the exact result 
based on numerically solving Eq. (|58|) . Plotted alongside 
this, although almost indiscernible, is the same solution 
expect using the nonexpanding spacetime kernel (dashed 
line), obtained by setting H —^ 0, a —^ constant in Eq. 
(|59|l . and the solution based on the adiabatic-Markovian 
approximation of Eq. (jSJ (dot-dashed line) for the same 
parameter set. As seen, the expanding and nonexpand- 
ing cases differ by very little and the adiabatic-Markovian 
approximation is in good agreement with the exact solu- 
tion. This result presented first in M confirms simplifying 
approximations claimed in a, 14, zO|, |23, 123 but up to 
now had not been numerically verified. 

More interestingly, the dotted line in Fig. ^is the so- 
lution that would be found by the conventional approach 
in which the nonlocal terms in E g. (I58II are ignored. The 
conventional approach 0, 0, [ifl 113 expects the inflaton 
to start oscillating, which is the precursor to entering 
various stages of pre/re-heating. However with account 
for dissipative effects, this never happens for our example 
in Fig. ^ since the inflaton remains overdamped till the 
end when it settles at its minima at 1^ = 0. Moreover, 
throughout inflation, and not just near the end, the infla- 
ton dissipates energy, which yields a radiation component 
of magnitude 



where a^ = N^h"^ /{%n). 



4H" 



(61) 
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IV. DENSITY PERTURBATIONS 

The difference between the cold and warm inflationary 
dynamics implies that the inflaton density perturbations 
in the two cases also have basic differences. In partic- 
ular in the warm inflation regime, since a thermalized 
radiation component is present with T > m^j^, inflaton 
fluctuations are dominantly thermal rather than quan- 
tum. There are two distinct regimes of warm inflation to 
note. One is the weak dissipative regime |24Ll25|. 

5ip^ - HT warm inflation (T < 3iJ), T > m^, (62) 

and the other is the strong dissipative regime [3], 

Sip"^ - VhTT warm inflation (T > 3H), T > m^. 

(63) 
For comparison, for cold inflation, where inflaton fluctu- 
ations are exclusively quantum l2fi| . 

Sifi"^ ^ H"^, cold inflation T < m^. (64) 

For both cold and warm inflation, density perturbations 
are obtained by the same expression, dp/p ~ Hdip/ip. 

The inherent difference in density perturbations in cold 
versus warm inflation provides a possible direction for 
distinguishing between these two inflation dynamics us- 
ing observation. In |23 an order of magnitude estimate 
of density perturbations during warm inflation was com- 
puted by matching the thermally produced fluctuations 
to gauge invariant parameters when the fluctuations cross 
the horizon (for other phenomenological treatments of 
warm inflation see l28l | ) . This work provided a clear state- 
ment of the consistency condition. Cold inflation has 
three parameters, related to the potential energy magni- 
tude Vb, slope en in Eq. (0) and t]h in Eq. 131, whereas 
there are four observable constraints {6h, Ag, Us, Ug). 
This implies a redundancy in the observations and allows 
for a consistency relation [23. This is usually expressed 
as a relationship between the tensor-to-scalar ratio and 
the slope of the tensor spectrum. Warm inflation has 
an extra parameter, the dissipation factor, which implies 
four constraints for four parameters. Hence we do not 
expect the consistency relation of standard inflation to 
hold in warm inflation 27] . Thus discriminating between 
warm and standard inflation requires measuring all four 
observables. The WMAP and upcoming Planck satellite 
missions should provide strong constraints on the scalar 
spectrum and having polarization detectors, it is hoped 
the tensor spectrum also will be measured. At the same 
level of approximation, nongaussian effects from warm 
inflation models were computed and found to be of the 
same order of magnitude as in the cold inflation case, and 
thus too small to be measured [Sfl ■ 

One interesting feauture about warm inflation dynam- 
ics is that it offers a solution to the ry-problem [sj]. In 
standard inflation models 0, Qj S] i where inflaton evolu- 
tion is damped by the term 3i?iyj, the slow-roll condition 
amounts to "qu ~ 1, which equivalently means the poten- 
tial can not have mass terms bigger than ~ H'^ip^. 



Since Supersymmetry suppresses quantum corrections, 
thus can preserve the tree level potential, it has been a 
central idea in realizing such flat inflationary potentials. 
Of course, since inflation requires a nonzero vacuum en- 
ergy density, inevitably SUSY must be broken during the 
inflation period, thus possibly ruining the desired degree 
of flatness in the potential. In particular, once supergrav- 
ity effects are included, it becomes very difficult for this 
symmetry to preserve flatness at the leevl oirjH < 1 ■ For 
F-term inflation, where the nonzero vacuum energy den- 
sity arises from terms in the superpotential, no symmetry 
prohibits the appearance of the Planck mass suppressed 



higher dimensional operators a„ t/?" /m"; |3j, |33|, |3J, |3 
For the large class of chaotic inflation type models |3], 
where inflation occurs with the inflaton field amplitude 
above TOpj, to control these higher dimensional opera- 
tors would require the fine-tuning of a infinite number 
of parameters or choosing only certain types of SUGRA 
corrections, such as the mininal Kahler potential. Even 
for models where inflation occurs for field amplitudes be- 
low the Planck scale, dimension six operator terms of the 
form Vcp"^ /m'ii ~ H^cp^ can emerge and ruin the desired 
fiatness. Both minimal and nonminimal Kahler poten- 
tials can lead to such terms J35:j ■ 

One possible solution to the yy-problem might be D- 
term infiation. In such models, the nonzero vacuum en- 
ergy arises from the supersymmetrization of the gauge 
kinetic energy. However a closer examination |34l l35j | re- 
veals that attaining the required degree of fiatness makes 
such models very restrictive. 

Up to now, attempts to solve the eta-problem have 
sought symmetries that can maintain this desired degree 
of flatness. One of the few that has proven successful 
is called the Heisenberg symmetry J33, although it is 
very restrictive. Another proposal has been a certain 
shift symmetry |35j . which is particularly interesting as 
it does not require SUSY. In common, all attempts so far 
have one foregone conclusion, that inflaton dynamics is 
only viable for tjh < 1. However, if the inflaton evolution 
happened to have a damping term larger than 3Hip, then 
clearly slow- roll can be satisfied for 77// > 1. Such a 
possibility is precisely what occurs in warm infiationary 
dynamics. 

To see this let us examine the warm inflation solution 
for the simple potential 



^-2< 



(65) 



In the cold inflation case, such a model requires an initial 
inflaton amplitude (0) = ip > nipi. Moreover, SUSY 
models that realize a potential like this inevitably lead 
to an eta-problem based on the reasons discussed above. 
Let us now treat this model in the warm inflation case. 
To focus on the essential points, our calculations here will 
be purely phenomenological, although they can be read- 
ily derived from a first principles quantum field theory 
calculation as done in Sect. IIIII We consider the case 
where the dissipative coefficient in Eq. (O is independent 



of both ip and T, T = constant. The inflaton initially 
is at a nonzero field amplitude (/9 ^ 0, thus supporting a 
vacuum energy. 

The background cosmology for models with constant 
T and monomial potentials has been solved exactly [g. 
From this we find Ne ~ HT/rn^. The radiation produc- 
tion is determined from the energy conservation equation 
(|14|l . During warm infiation p^ ~ 6, 25], so that Eq. 
(|14|l reduces to Eq. H61|l . Identifying pr ~ T'^ permits 
determination of the temperature during warm infiation. 
Finally, once T is determined, Eq. (|63|l allows determi- 
nation of density perturbations. 

Combining these expressions, for model Eq. (|65|l with 
T = const, in Eq. Q, gives 



iVe ~ 2\/2- 
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(66) 



(67) 



(68) 



Imposing observational constraints N^, = 60 and 6p/ p ■ 
10^^, leads to the results 



H niff, 



^ « 5.3 X 10«^ 



nipi 



■nipi 



(69) 



(70) 



T « 4x 10 ^nipi, and T « 10 ttt.^, with the ratio niip/mpi 

free to set. For tti^/ttt-p/ ~ 10~^, it means rju > 1 
and if < mpi. Thus we see for sufficiently small infia- 

ton mass, m,^ ~ 10^"GeV, there is no eta-problem, since 
m^ ^ H and tp < nipi. Since this warm inflation so- 
lution works for rjH 3> 1, SUSY models realizing simple 
monomial potentials like Eq. H65I) do not require any 
special symmetric as is the case for cold inflation models. 
The "eta" and large (^-amplitude problems simply cor- 
rect themselves once interactions already present in the 
models are properly treated. 



V. HYBRID SUSY MODEL 

Following the analysis of 12i|, Let us apply the results 
of the previous two sections to the SUSY hybrid model 
matter held A A coupled to it 



W = k5($i$2 - M^) + .g$2AA . 



(71) 



In this model, the inflaton is identifled with the bosonic 
part of S, 4>s. The above is a toy model representing 



an example of how the basic hybrid model, first term 
on the RHS, is embedded within a more complete par- 
ticle physics model, in this case through the A fields. 
We will show that in the above model both cold and 
warm infiation exist and we will determine the parame- 
ter regime for them. This will then explicitly verify the 
conclusions from the recent papers on dissipation [3, |31 j 
that showed both types of infiationary dynamics could 
exist. In both infiationary regimes, we will calculate the 
scalar spectral index ns — 1, and its running dns/dhik. 
With this information, we will then identify the qualita- 
tive and quantitative differences arising from the warm 
versus cold regimes. 

In this model the heaviest field with mass mj^ can de- 
cay into the massless fermionic partners of A and A, with 
decay rate: 



167r 



(72) 



We have F i 



J , and this can be much larger than the 



Hubble rate during inflation: 
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(73) 



where x^ = (f>s/{V^p)- Having F^/iy > 1, all the way 
up to the end of infiation, only requires g > 0.16 for k < 
0.001 ( g > 0.01 for k = 0.5). This allows us to apply the 
adiabatic-Markovian limit in the effective EOM for the 
inflaton background field, Eq. (0, with the dissipative 
coefficient given now by: 
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and the ratio to the Hubble rate is given by: 
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(74) 



(75) 



which behaves like T/(3iJ) ex xn esc (j)s, and so decreases 
during inflation. That is, the evolution of the inflaton 
field may change from being dominated by the friction 
term T to be dominated by the Hubble rate H. Whether 
the transition between these two regimes happens before 
or after 60 e-folds will depend on the value of the parame- 
ters of the model like k and g. The amount of "radiation" 
obtained through the dissipative term, is given by: 
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(76) 



which even when T < H could give rise to a thermal 
bath with T > H. In particular, we can have: 
(a) T > SH, and T > i? ((^s ^ -V^^/T): 
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Note one might expect the presence of this thermal bath 
may induce thermal corrections to T and T4// but as 
shown in 36], these corrections are negligible. 

The values of the couplings k and g for which we could 
have cold or warm inflation, and strong or weak dissipa- 
tive dynamics, are plotted in Fig. ^. In order to get the 
different regions, we have proceeded as follow: for each 
pair of values in the plane n — g, the value of the inflaton 
field at the end of inflation is determined. This is done in 
the cold and weak dissipative regimes by the condition^ 
r^x = 1, Eq.Q. In the strong dissipative regime inflation 
can end either with 7]x = 1 or it may also happen that 
most of the vacuum energy is already transferred into ra- 
diation during inflation, and then inflation will end when 
pn ~ K^n'^ instead. In this case, whichever occurs first 
fixes the value of the inflaton field at the end of inflation. 
The value of the inflaton field at 60 e-folds of inflation is 
then obtained from 



iV, 



AV 



;(l + r)d0. 



(79) 



This in turn fixes the value of the dissipative coefficient 
T, Eq. (|74|l . the temperature of thermal bath, Eq. (|76|l . 
and therefore the amplitude of the spectrum P-ji. The 
COBE normalization is then used to fix the value of the 
scale /z. In order to match the expressions for the spec- 
trum across the different regimes, we have used a simple 
expression with : 



pl/2 

n 




(80) 



We can see in Fig. that the strong dissipative 
regime T > SH requires large values of the couplings, 
K r^ g ^ 0(1); for values k c^ g "^ 0.1 we are in the weak 
dissipative regime; and for values k ~ g ~ 0.01 we recover 
the cold inflationary scenario. Typically, for a fixed value 
of the scale p. the amplitude of the spectrum in the strong 
dissipative regime would be larger than the one gener- 
ated at zero T . The COBE normalization implies then a 
smaller value of the inflationary scale fx. For example, for 
K = g = 1 we have /j, ~ 10"'^'^ GeV, whereas pushing the 
coupling toward its perturbative limit, k = g = VItt we 
would get /x ~ 2 X 10^^° GeV. On the other hand, going 
from the cold to the weak dissipative regime, the value 
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FIG. 2: Regions of cold {T < H), and warm (T > H) SUSY 
hybrid inflation in the k, — g plane. The warm inflation re- 
gion is divided into the weak dissipative regime with T > H 
and Ts < 3// (lighter shaded region) , and the strong dissipa- 
tive regime with T5 > 3// (darker shaded region). Included 
are also the contour plots of constant ^, and the adiabatic- 
Markovian limit r+ = H . The black region on the right of 
the plot is excluded because (^5 > mp. In addition, when 
SUGRA corrections are taking into account, values to the left 
and down the wide dot-dashed line are excluded. 



of /i only varies by a factor of 2 or 3, and it is still in the 
range of the GUT scale 0{lQ^^) GcV. 

In Fig. © we have compared the prediction for the 
spectral index of the scalar spectrum of perturbations in 
both the CHI scenario, and warm hybrid inflation (WHI). 
From the warm inflation scenario we can always recover 
the CHI prediction by taking g <^\. In standard SUSY 
GUT hybrid inflation, for small values of the coupling 
K the spectrum is practically scale invariant, it reaches a 
minimum around k ~ 0.01, and then rises due to SUGRA 
corrections up to positive values, which are disfavoured 
by WMAP results. But in the weak and the strong dissi- 
pative regime, due to the different origin of the spectrum, 
we get that the spectral index is still below 1 even for 
values of the couphng n > 0.01. Specially in the strong 
dissipative regime, where the dynamic is such that the 
inflaton field is well below the Planck scale and SUGRA 
corrections are negligible. In that regime the departure 
from scale invariance is within the observational value, 
with ns-l~ -0.022. 

As a well-motivated example |l2| . which combines in- 
flation with leptogenesis and light neutrino masses [33, 
l3g |. the inflaton can decay into right handed (s)neutrino 
fields VRi (i= family index). The decay proceed through 
the non-renormalizable coupling ^i^ii'mi'B.i, with decay 
rate, 



^ The value of ryx becomes larger than 1 before the other two slow- 
roll parameters. 
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FIG. 3: Spectral index for cold SUSY hybrid inflation (solid 
line, CHI), and warm inflation (gray region, WHI). The weak 
dissipative regime {T > H but T s < 3//) is given by 
the darker gray region (triangle down);the strong dissipative 
regime (Tg > 3-ff) is given by the light gray region (triangle 
up). 



flation picture. An important point that has been at- 
tempted to be converyed is that interactions in a infla- 
tion model Lagrangian can have significant effects during 
the inflation period. This point has led to the key result 
that many typical models of inflation, which had been as- 
sumed to yield just cold inflation dynamics, in fact have 
regimes of warm inflation. The realization of this fact 
leads to many interesting and new questions. From the 
one direction is model building. In the warm inflation 
case, there are some new features that do not exist for 
cold inflation. For one the ry-problem can be elminated 
simply by the dissipative dynamics and does not require 
any further constaints on the particle phyiscs model, in 
particular the Kahler potential. This means greater flex- 
ibility to the high energy properties of the models. An- 
other new feature is that for monomial potentials, ob- 
servationally consistent warm inflation occurs for field 
amplitudes below the Planck scale, in constrast to cold 
inflation, where it is well known chaotic inflation requires 
if > TTipi. This difference implies from the perspective of 
effective theories, monomial potentials are accpetable in 
the strong dissipative warm inflation regime. 



where M.^ is the RH (s)neutrino mass. In the CHI sce- 



nario, with fi ~ 0(10^^) GeV, n ~ 10"^, and mg ~ 10 
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GeV, the gravitino constraint Tjih < 10^ GeV trans- 
lates roughly into M, ~ 10" V ~ O{10^^)GeV. Those 
values are also consistent with baryogenesis and light 
neutrino masses 38]. This kind of scenario is also vi- 
able in the warm inflationary regime. Being consis- 
tent with the observed baryon asymmetry and the atmo- 
spheric neutrino oscillations does not directly constraint 
the value of k but the value of ms ~ 10"'^^ GeV. In the 
warm inflationary regime the value of the scale ^ re- 
quired for successful inflation reduces as we moved into 
the strong dissipative regime, ms is of the order of lO^'^ 
GeV for k ~ 0(1), and the gravitino constraint gives 
now Mi ~ 10-3/1 - 0(101") GeV. Therefore, a model 
of warm inflation and leptogenesis without the need of 
small couplings would be viable and compatible with ob- 
servations, in the strong dissipative regime. 



Another interesting direction is determining tests that 
could discrimnate warm versus cold inflation from obser- 
vational data. The nature of the density perturbations 
in the two dynamics is different, so one should expect to 
find different observational signatures in the two cases. 
For one the tensor-to-scalar rations differ for the two 
cases. Also, although we have not explored this direc- 
tion in great detail in this talk, more careful analysis of 
density perturbations done by evolving the complete set 
of cosmological perturbation equations was done in 221 ■ 
This analysis showed that dissipative effects can produce 
a rich variety of scalar spectra ranging between red and 
blue and can also cause oscillations in the scalar power 
spectrum. Further understanding of all these possibilities 
is needed before firm criteria can be set of for discrimi- 
nating between warm and cold inflation behavior in the 
CMB data. 



VI. CONCLUSION 

In this talk we have reviewed the basic ideas of warm 
inflation and compared them to the standard or cold in- 
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